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Abstract 

The total cross section of the pd — s- pdrj reaction has been measured at 
two energies near threshold by detecting the final proton and deuteron in a 
magneti spectrometer. The values are somewhat larger than expected on the 
basis of two simple theoretical estimates. 
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The production of rj mesons in the pd ^ ^He r] reaction near threshold is remark- 
able for both its strength and its energy dependence |I|, Q . The threshold amplitude 
is of a similar size to that for pion production, despite the much larger momentum 
transfers associated with i] formation. Although the angular distribution remains 
isotropic, suggesting 5'-wave production, the square of the amplitude falls by a factor 
of three over a 5 MeV change in the cm. excess energy Q. This has been taken as 
indicative of a nearby quasi-bound state of the t] ^He system, arising through strong 
T) multiple scatterings from the three nucleons in the recoiling nucleus 0. 

In order to transfer such large momenta, Kilian and Nann suggested two-step 
processes involving intermediate pions. They showed that the threshold kinematics 
for pd ^ ^He 7] were in a sense magic. The momentum of the t] produced in the 
reaction is very similar to that obtained from the sequential physical processes of 
pp diT^ followed by n^n prj, when there is no relative momentum between the 
final pd pair and all Fermi momenta are neglected. In such cases the final proton 
and deuteron are likely to stick to form the observed ^He nucleus. The classical 
estimate of the enhancement due to the magic kinematics is broadly confirmed by 
quantum mechanical calculations , which reproduce the size of the near-threshold 
cross section to within about a factor of two. 

The same two-step model should also be capable of explaining events where the 
final proton and deuteron emerge freely in the pd pdrj reaction and the aim of the 
present investigation was to undertake a first exploration of this cross section near 
the threshold beam energy of Tp = 901.2 MeV. Unfortunately, due to the closure of 
the laboratory, data could only be obtained at two beam energies. 

The experiment was carried out at the Laboratoire National SATURNE (LNS), 
using the large acceptance magnetic spectrometer SPESIII, which was well adapted 
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for the study of meson production in three-body final states near threshold through 
the detection of two charged particles. The experimental conditions concerning the 
beam monitoring, particle detection and identification, were rather similar to those 
of previous studies of meson production in the pp ^ pp X , where the meson X was 
identified by the missing mass method |p. A liquid deuterium target of 207 mg/cm^ 
thickness was employed and, in order to improve the missing mass resolution and 
the signal-to-background ratio, the opening of the vertical collimators of SPESIII 
was reduced to ±40 mr. 

One special feature of the pd pdr] reaction near threshold is the rather low 
momentum, around 400-500 MeV/c, of the outgoing protons. This is to be compared 
with the standard 600-1400 MeV/c momentum range of the SPESIII spectrometer. 
The momentum of the recoiling deuterons is about 900 MeV/c and, in order to detect 
both particles simultaneously, the magnetic field was tuned down to accept momenta 
from about 360 to 960 MeV/c. Under normal SPESIII working conditions, the 
values of the particle momenta were obtained by using well established polynomial 
relations, taking the coordinates of the trajectories near the focal surface as input. 
The properties of SPESIII were not extensively studied with reduced fields and, in 
the present experiment, we used the polynomial parametrisation with the momenta 
of the particles scaled according to the ratio of the actual to the standard mean 
field, (2.03 Tesla)/(3.07 Tesla). This procedure essentially assumes that the field 
was reduced uniformly. A similar method was applied in the simulations, applying 
the same ratio to the momenta when tracking the particles. Such simulations are 
important for generating the expected missing mass peak of the pd —>■ pdrj reaction 
as well as the background spectrum of the pd ^ pd2TT reaction. 

Two-dimensional experimental and simulated scatter plots of the emerging pro- 
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ton and deuteron momenta are shown in Fig. 1. Superimposed upon a fairly uniform 
background, due the pd ^ pd27T reaction, there is a darker elhpse inside which the 
pd —>■ pdrj events are confined. The experimental missing mass spectra of the 
pd pdX reaction are shown in Figs. 2al and 2a2. Clear 77 peaks are observed 
near the upper edges of phase space for the two nominal proton beam energies of 905 
and 909 MeV. Simulated background spectra of the pd ^ pd27r reaction are shown 
in Figs. 2bl and 2b2 and simulated peaks of the pd ^ pdrj reaction are represented 
by the solid lines in Figs. 2cl and 2c2. To evaluate the number oipd —>■ pdrj events, 
the two simulated spectra were combined so as to fit the experimental data. After 
subtracting the simulated background from the experimental results, the remaining 
events (points with error bars) in Figs. 2cl and 2c2 show good agreement with the 
simulated pd ^ pdr] spectra (sohd line). 

Taking the mass of the rj meson to be 547.30 MeV/c^ |^, the best fits were 
obtained by assuming incident proton energies which were 1-2 MeV lower than the 
nominal values derived from the Saturne machine parameters. The fits also sug- 
gested adjusting the mean field ratio to a value slightly below that of the initial 
2.03/3.07 ansatz. However, due to the uncertainties in the effective field strength 
and the particle tracking, no definitive accurate values of the beam energies could be 
deduced from the fitting procedure. But, since the shift indicated here is very similar 
to the mean difference AT = T„o^i„^i — Treasured = 1.1 ± 1.0 MeV obtained at Saturne 
from other meson-production reactions near threshold ^, ^ , we adopt this energy 
correction AT, to be subtracted from the nominal values to obtain the 'true' ones. 



The average values of the excess energies Q = y {mp + rriciY + 2mdTp—mp—md—m^, 
where the rrii are particle masses, were determined using the corrected proton ener- 
gies and taking into account energy losses in the target. These energies were used in 
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the simulations required to evaluate the SPESIII acceptances, which were estimated 
assuming phase-space distributions of final particles. The acceptance decreases very 
quickly above threshold through one of the final particles falling outside the solid 
angle of SPESIII. Nevertheless, from the resulting angular dependence of the ac- 
ceptance shown in Fig. 3 as a function of the cosine of the t] emission angle in the 
centre-of-mass system, it can be seen that all the regions of phase space were cov- 
ered. The resulting overall acceptances of 24% and 7%, evaluated at the mid-target 
energies of Tp = 903.7 and 907.7 MeV respectively, lead to the cross sections given 
in Table 1. The first error on the cross sections given in the table includes the statis- 
tical error (10% and 18% respectively) and a 13% systematic error on the absolute 
normalisation. The ±0.6 MeV uncertainty in Q gives rise to the additional quoted 
errors through the rapid acceptance variation. However, the latter errors affect little 
the comparison with theory since, if both values of Q are increased by 0.6 MeV, the 
experimental points move largely in the directions given by the theoretical curves. 
Estimates have been made of the pd — pdrj total cross section near thresh- 



old in the quantum two-step model of Ref. |T0|, though neglecting all final state 
interactions. The energy variation of 

axipd^pdri) = 1.2 nh (1) 

is compatible with that of our data shown in Fig. 4. The predicted values of 2.7 and 
17 nb, at Q = 1.5 and 3.8 MeV respectively, are only about a factor of two smaller 
than our results and this discrepancy could be due to the neglect of the strong final 
state interaction between the proton and deuteron. 

There are two possible 5- wave proton-deuteron final states, corresponding to spin 
I and |. The low energy spin-quartet scattering wave functions show little structure 



at short distances, whereas the spin-doublet bear some similarity to the shape of the 
pd distribution at short distances inside the bound ^He nucleus ||TT|. 

The connection between the bound and scattering pd wave functions can be 
exploited to estimate the production of r] mesons in the pd ^ pdr] reaction, with 
the final pd system in the spin-| state, in terms of the cross section for pd ^ ^Herj. 
Provided that the structure of the deuteron is neglected, the relative normalisation 
of the bound and scattering wave functions at short distances is fixed purely by the 
proton-deuteron binding energy, e ^ 5.5 MeV, in the ^He nucleus. If the meson 
production operator is also of short range, the production in two and three-body 
final states should be related through |12| 

aripd ^ pdr]) = \ [j) ^ (l + \/l + ) ' x aripd ^ ^He r]) . (2) 

This approach reproduces about | of the pd pdii^ total cross section [jl3 
in terms of that ioi pd ^ '^He7r° and the residue could be due to spin-| final 
states. 

The very precise near-threshold pd ^ ^He rj total cross section data may be 
parametrised as 

,,(p,^3He„)^g)^j^^^J^^,b. (3) 

where the rj and proton cm. momenta and Pp are measured in fm~^. This 
parametrisation is shown in Fig. 4. 

The predictions of Eq. (2) for the pd ^ pdrj total cross sections are about a factor 
of three lower than our experimental results shown in Fig. 4. This may be due to the 
short-range assumption for the meson-production operator made in deriving Eq. (2). 
In the two-step model of ref. 0, the momentum transfer is provided through having 
a secondary interaction with an intermediate pion. Since high Fermi momenta are 
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then not required, this means that the final p d system is not necessarily produced 
at short distances. However, at larger distances the scattering wave functions are 
generally bigger than bound state ones, which must die off exponentially. It would 
therefore be very desirable to have a microscopic two-step model calculation of the 
type of ref. but with the proton-deuteron final state interaction included. 

We have made the first measurements of the pd ^ pdr] reaction near threshold 
and obtained cross sections about a factor of 2-3 higher than those of two simple 
theoretical approaches. Since, for our data, Q is less than the ^He binding energy e, 
the difference in the energy dependence of the two present models comes principally 
from the striking behaviour of the pd — > ^He rj cross section. More detailed exper- 
iments, with a much better resolution in Q, are required to see if there is in fact a 
strong 7] final state interaction in the pdi] system to match that in ^B.er]. 

We wish to thank the Saturne accelerator crew and support staff for provid- 
ing us with working conditions which led to the present results. Discussions with 



U. Tengblad regarding Ref. |10| were very useful. 
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Table 1: Measured total cross sections for the pd ^ pdrj reaction at a mid-target 
centre-of-mass excess energy Q. Although there are ±0.6 MeV uncertainties in the 
values of Q, the relative value is correct to ±0.1 MeV. The first quoted error in 

the cross sections includes statistical and normalisation uncertainties; the second is 
that induced by the uncertainty in Q. For the first point the energy uncertainties arc 
negligible for a/Q, whereas for the second point it is a/Q"^ which is largely unaffected 
by Ag. 



Q 

(MeV) 


aripd — > pdrj) 
(nb) 


1.5±0.6 


6.1±1.0±|8 


3.8±0.6 


40±9±ig 
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Figure 1: Two-dimensional scatter plot of the momenta of the proton and deuteron 
arising from the reaction pd ^ pdX a. little above the rj threshold. The experimental 
data, shown in the upper figure, were measured at a nominal beam energy Tp — 
909 MeV but, as discussed in the text, this is expected to be 1-2 McV above the 
true value. The simulation in the lower figure was carried out at 907 MeV. Events 
from the pd ^ pd2'ri reaction populate the whole plot but in the top left corner one 
can see an ellipse containing extra events corresponding to rj production. 
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Figure 2: Missing mass spectra of the pd ^ pdX reaction at nominal beam energies 
of (1) 905 MeV, and (2) 909 MeV. The observed spectra in (al) and (a2) are to be 
compared with background spectra prf ^ pd2'K m (bl) and (b2) estimated by Monte 
Carlo simulation. The shapes and widths of the rj peaks, found by subtraction and 
shown in (cl) and (c2), are in good agreement with the Monte Carlo predictions (solid 
lines) assuming true beam energies 1-2 MeV below the nominal values suggested by 
the macroscopic parameters of Saturne. 
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Figure 3: Acceptance of the SPESIII spectrometer for the reaction pd ^ pdrj aX 
mid-target beam energies of 903.7 and 907.7 MeV as a function of the cosine of the 
77 production angle in the centre-of-mass system. The curves have been calculated 
assuming phase-space distributions. 
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Figure 4: Measured total cross sections for the pd ^ pdrj reaction close to threshold. 
The additional errors arising from the influence of the uncertainty in the absolute 
value of Q given in Table 1 are not indicated since, due to the Q — a error correlation, 
they have much smaller components perpendicular to the theoretical predictions. 
The solid curve is the prediction for the spin-| final state in the final-state-interaction 
model of eq. (2) [|l^]. The dashed curve has been estimated in the two-step model 
of ref. []TU|, where final state interactions have been neglected. The pd ^Het] 
data of ref. 0] and the parametrisation of Eq. (3) are shown by the open circles and 
dot-dashed curve respectively. 
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